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ABSTRACT 

We present near-infrared (NIR.) time-series spectroscopy, as well as complementary ultraviolet (UV), optical, and NIR. 
data, of the Type la supernova (SN la) iPTF13ebh, which was discovered within two days from the estimated time of 
explosion. The first NIR spectrum was taken merely 2.3 days after explosion and may be the earliest NIR spectrum 
yet obtained of a SN la. The most striking features in the spectrum are several NIR Cl lines, and the Cl A1.0693 //m 
line is the strongest ever observed in a SN la. Interestingly, no strong optical C II counterparts were found, even though 
the optical spectroscopic time series began early and is densely-cadenced. Except at the very early epochs, within a 
few days from the time of explosion, we show that the strong NIR CI compared to the weaker optical CII appears to 
be general in SNe la. iPTF13ebh is a fast decliner with Amis(B)= 1.79 ± 0.01, and its absolute magnitude obeys the 
linear part of the width-luminosity relation. It is therefore categorized as a “transitional” event, on the fast-declining end 
of normal SNe la as opposed to subluminous/91bg-like objects. iPTF13ebh shows NIR spectroscopic properties that 
are distinct from both the normal and subluminous/91bg-like classes, bridging the observed characteristics of the two 
classes. These NIR observations suggest composition and density of the inner core similar to that of 91bg-like events, 
and a deep reaching carbon burning layer not observed in slower declining SNe la. There is also a substantial difference 
between the explosion times inferred from the early-time light curve and the velocity evolution of the Sill A0.6355 /im 
line, implying a long dark phase of ~ 4 days. 

Key words, infrared: general; supernovae: general; supernovae: individual: iPTF13ebh 


1. Introduction 


Type la supernovae (SNe la), with empirically calibrated 
luminosities, provide a direct measure on the expansion 
history of the univers e and led to the discovery of the 
accel erated expansion (Ri ess et al.l 119981 : iPerlmutter et all 
19991) . The goal of re aching the 1 — 2% distance precision 
(IAlbrecht et al.l[2006D or o w ss 2er m dGoobar fe Leibundgutl 
12011 ) for the next-generation SN la experiments critically 
depends on minimizing the systematic errors in the distance 
determinations. One of the most important is the potential 
evolution of the empirical calibration with redshift. Given 
the large range of look-back time considered, such an evolu¬ 
tion could be caused by changes in th e mean host m etallic - 
ity and stellar mass with redshift (e.g. JSullivan et aHl2010h . 
A deeper understanding of the physics of SNe la could help 
to estimate this systematic error and perhaps even mitigate 
its effects. 


* This paper includes data gathered with the 6.5-meter Mag¬ 
ellan Telescopes located at Las Campanas Observatory, Chile. 


The general consensus for the origin of a SN la is the 
thermonuclear explosi on of a carbon-oxygen white dwarf 
(IHovle fc Fowleiill96 0) . Since oxygen is also produced from 
carbon burning, carbon provides the most direct probe of 
the primordial material from the progenitor. The quantity, 
distribution, and incidence of unburned carbon in SNe la 
provide important constraints for explosion models. Turbu¬ 
lent deflagration models predict that a large amount of un¬ 
processed carbon should be left over dGamezo et al.l [2003 : 
iRonke et al.l [20071. Three-dimensional simulations of pure 
deflagrations predict a significant amount of unburned ma¬ 
terial in tire inner ejecta, in conflict with observations (e.g., 
iKozma et al.ll2005t l . In one-dimensional simulations, a tran¬ 
sition to a deton ation is expected t o result in complete 
carbon burning (IHoflich et al.l l2002t iMarion et alJ 120061 : 
iKasen et al.l l2009h . although there are exceptions. Three- 
dimensional effects in delayed detonation mo dels can create 
pocke ts of unburned material in the ejecta dGamezo et al.l 
I2004D . In the case of puls ating delayed detonation models 
(PDD; IHoflich et al.lll996h . the initial pulsation allows the 
outer ejecta to become loosely bound. The low density of 
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this material subsequently quenches the nuclear burning. 
This leaves more unburn ed carbon than the standard de¬ 
layed detonation models dDessart et al.|[2014l) . 

The identification of the weak absorption feature near 
0.63 pm in the normal SN la SN 199 8aq as C II AO.6580 pm 
was fi rst suggested bv lBranch et al.1 (1 20031) . iGaravini et al.l 
(|2004h then noted possible contri butions of C II and C III in 
SN 1999aa. [Thomas et al.l (l2007f ) presented convincing de¬ 
tections of the Cn AA0.4267, 0.6580, and 0.7234 lines in 
the early optical spectra of SN 2006D. Cn A0.6580 pm 
has also been detected in the luminous and slowly ex¬ 
panding super- Chandrasekhar candid a tes: strong C II lines 
in SN 2009dc dSilverman et al.l 12011 ;: Taubenbererer et al 


2011), and marginal de tections in SN 2003 fg ([Howell et al 


2006) and SN 2007if ( Scalzo et al.l l2010h . In subsequent 
studies with larger optical spectroscopic samples, 20 — 
30% of the pre -maximum spectra were found to show 
Cll signatures ( iThomas et al.l 20 lit Parrent et ahl 120111 : 


120061 : iThomas et al.ll2007r Tanaka et al.1 20081) . SNe la that 


sion physics (e. 

2.. IWheeler et al. 19981: Hoflich et al. 2004; 

Motohara et al. 

20061: iFriesen et al.l 2014 iDiamond et al.1 


file. The NIR CI line has the same velocity as the optical 
C II A0.6580 pm line and increases in strength toward maxi¬ 
mum light, in contrast to the fast-fading nature of the opti¬ 
cal C II line. iHsiao et all (120131 ) suggested that this delayed 
onset of the NIR CI line is an ionization effect and pointed 
out its potential to secure more representative properties of 
unburned material in SNe la. 

iPTF13ebh is a SN la discovered at an exceptionally 
young age by the intermediate Palomar Transient Fac¬ 
tory (iPTF). The first NIR spectrum was taken merely 
2.3 days past the explosion, allowing the examination of 
material in the outermost ejecta. Strong NIR Cl A1.0693 
pm is present in the first two NIR spectra. To our knowl¬ 
edge, this marks the fourth detectio n of NIR Cl so fa r 


in a S N la, along with S Ne 1999by (Hoflich et al 


201 lfe (IHsiao et al.ll2013l) . and 2014J ( Marion et al.l 


2002 ). 


2015). 


IFolatelli et al.|[2012l : Silverman fe Filippenkoll2012fl . Several 
of these studies noted, however, that this fraction represents 
a lower limit, as noise, line overlap, and the phase at which 
the SNe la were observed could affect the C n detection 
(e.g.. iBranch et aLll2007l : IFolatelli et al.1120121) . 

The velocities of the detected C II lines are gener¬ 
ally low. There have also been hints that SNe la with 
detected C II hav e preferentially bluer colors and nar¬ 
rower light curves (|Thom as_eta |2014|jR)latejIi et_al.l l2012l : 


[Silverman fe Filippenkol 120121 Maguire et al.l 2014T ) . but 

this preference was not clear in the examinations of othe r 
data sets (e . g.. |Parrent et al.l 120111 : iBlondin et al.l 120121) . 

I Milne et al.l (120131) also noted a strong preference for blue 
ultraviolet (UV) colors in SNe la with carbon detections. 
There appears to be a consensus that the mass fraction 
of carbon is low (e.g., Branch_et_aL 20031: iMarion et all 


display signatures of carbon tend to have lower Sill ve¬ 
locities, while the objects without carbon span the en¬ 
tire range of Sill velocities (e.g., IFolatelli et al.l 120121 : 
ISilverman fe Filippenkol I2012D . This phenomenon may be 
an observational bias as the CII A0.6580 pm at high veloc¬ 
ities is shifted into the prominent Sill A0.6355 pm line. 
It cou ld also be a physical bias. The PDD models of 
iDessart et al.l (|2014l) predict both strong carbon lines and 
relatively narrow Sill A0.6355 pm absorption, while their 
standard delayed-detonation models predict no carbon and 
broad Si II lines. Several characteristics of unburned mate¬ 
rial established by the optical CII lines may be uncertain 
because of the observational challenges posed by capturing 
an unbiased sample of these features. 

Near-infrared (NIR) spectroscopy of SNe la has been 


In fact, the Cl feature of iPTF13ebh is the strongest ever 
observed. As the supernova evolved, it became clear that 
the time evolution of the NIR Cl feature was very differ¬ 
ent from that of any other previous detections, and that 
iPTF13ebh had fast-declining light curves. Examination of 
the optical spectra showed no apparent T ill feature that 
woul d qualify it as a 91bg-like object dFilippenko et al.1 
fT99l . The NIR primary maxima also occurred before 
the B maximum. This places iPTF13ebh in the class 
of “transitio nal” objects, whose notable members include 
SNe 1986G (iPhillips et al.lll9 87), 2003gs (iKrisc iunas et all 
120091) . 2 004eo J Pastorel lo et al.l 2007 ^ Mazzali_et_ak[ 20081) . 
2009an |Sahu^_aL 2013 ), 201 liv ( Fblev et al.l 2012 ). and 
2012ht (|Yamanaka et al.1 20141) . In this paper, we focus on 
the examination of the NIR spectroscopic properties of 
iPTF13ebh, and highlight how they differ from those of 
normal and 91bg-like objects. These differences could shed 
light on the connections between normal and 91bg-like ob¬ 
jects. 


2. Observations and Reduction 


iPTF13ebh was discovered 
i mage taken with the 


an 
(P48 


carbon features. Examining the Cl lines in NIR spectra 
of three normal SNe la, they concluded that the abun¬ 
dance of unprocessed material is low. The early discovery 
of the nearby SN 201 lfe offered an unprecedented oppor¬ 
tunity to study the NIR carbon features in detail. With 
the aid of high signal-to-noise NIR spectr a and the au¬ 
toma t ed spectrum synthe sis code SYNAPPS (IThomas et al.1 
120111) , IHsiao et al.l (120131) inferred the presence of the CI 
A1.0693 pm line in SN 201 lfe from the flattening of the 
emission component of the Mgn A1.0092 pm P-Cygni pro¬ 


in NGC 890 by iPTF in 

_ _ Palomar 48-inch telescope 

Rahmcr et al. 120081) on 2013 November 13.15 UT 
dGao et al.l 12013(7 ! Within hours of the discovery, the 
Carnegie Supernova Project (CSP) obtained a NIR spec- 
trum with t he F olded-port Infrared Echellette (FIRE; 
ISimcoe et al.l 120131 ) mounted on the Magellan Baade Tele¬ 
scope on 2013 November 14.16 UT. Using a quick reduc¬ 
tion pipeline, based on the firehose packaged, a spec¬ 
trum was produced from the first two frames within 10 
minutes of acquiring the target. The supernova was clas¬ 
sified as a young SNIa with several str ong features pro - 
duced by intermediate-mass elements (ICao et al.1 l2013f ). 
More frames were subsequently obtained to reach the de¬ 
sired signal-to-noise ratio. An optical spectrum taken on 
2013 November 15.83 UT later confirmed the NIR clas¬ 
sification dMaguire et al.l [20131) . Reprocessing of the P48 
images from preceding nights showed that the supernova 
was visible on an image taken on 2013 November 12.21, at 
20.9 ± 0.2 mag, and not detected in an image taken 2013 
November 11.25 UT down to 21.7 mag. Intense follow up 
of iPTF13ebh was conducted in the optical and NIR as 
part of the CSP, in collaboration with iPTF and the CfA 
Supernova Group. 


Available at http://web.mit.edu/~rsimcoe/www/FIRE/. 
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Hsiao et al.: NIR carbon of SN la iPTF13ebh 


Table 3. Journal of du Pont RetroCam photometric observa¬ 
tions 


iPTF13ebh 



1 arcmin 



Fig. 1. A Swope r-band image of iPTF13ebh in NGC 890 taken 
near maximum light. The cross hair marks the location of the 
supernova. The compass and the size of the field are also noted. 


Table 1. Journal of P48 photometric observations 


MJD 

r 

56608.21 

20.94 (0.21) 

56609.15 

18.89 (0.06) 

56623.09 

15.14 (0.01) 

56624.08 

15.15 (0.01) 

56625.26 

15.15 (0.01) 

56626.14 

15.16 (0.01) 

56648.09 

16.60 (0.01) 


Notes. The MJD column lists the Modified Julian Date of each 
observation. The P48 photometry is host galaxy subtracted. 


2.1. Photometric observations 


The discovery and a small number of follow up images were 
taken with the P48 in the r band. Host galaxy subtrac¬ 
tions were performed using a deep reference of NGC 890 
constructed from images taken before the explosion. The 
point-spread function (PSF) was determined and matched 
prior to each subtraction, and was subsequently used for 
PSF photometry on the subtracted frames. The P48 pho¬ 
tometry is presented in Table [T] 

Swope follow up in uBVgri began on the same night 
that the FIRE classification spectrum was obtained. All 
images were taken with the newly commissioned e2v CCD 
imager which is more efficient overall and specifically more 
sensitive in the blue than the previous CCD employed by 
the CSP. The pixel size of 0'.'435 remains the same. The 
2x2 CCD array yields a field of view of approximately 
30' on each side. The supernova and standard fields are 
normally placed at the center of one pre-selected quad¬ 
rant. The bandpass functions with the new e2v imager have 
been characteri zed with a spec t romet er, in the same manner 
as des cribed in lRheault et all (|2010l ) and Stritzinger et~ahl 
(1201 ll) . The results are presented in lRheault et ali i 20141) . 
A Swope r-band image taken near maximum light is shown 


MJD Y J H 


56611.14 16.913 (0.016) 

56613.12 16.223 (0.011) 

56617.11 15.338 (0.009) 

56622.10 15.073 (0.009) 

56623.14 15.139 (0.009) 

56624.07 15.225 (0.009) 

56638.05 15.186 (0.009) 

56648.03 15.267 (0.009) 

56650.04 15.422 (0.009) 

56653.03 15.605 (0.010) 

56654.04 15.655 (0.011) 


16.643 (0.022) 
16.156 (0.014) 

15.059 (0.009) 
15.058 (0.009) 
15.196 (0.009) 
15.996 (0.012) 
16.317 (0.014) 
16.552 (0.019) 
16.884 (0.023) 
16.956 (0.028) 


16.627 (0.028) 
16.063 (0.016) 
15.315 (0.011) 

15.206 (0.011) 

15.207 (0.013) 
15.303 (0.011) 
15.291 (0.012) 
15.824 (0.022) 
16.063 (0.020) 

16.277 (0.035) 


Notes. The MJD column lists the Modified Julian Date of each 
observation. The RetroCam photometry is without host subtrac¬ 
tion. 


in Fig. [I] The reduction of Swop e images was done as 
described in IContreras et all ( 2010 ). PSF photometry was 
performed with re spect to a local seque nce of standard stars 
calibrated to the lLandoltl (|l992f) and ISmith et al.1 (|2002fl 
standard fields. The standard fields were observed on the 
same nights as the supernova observations over the course 
of 20 photometric nights. The Swope photometry is tabu¬ 
lated in Table [2j 

The Swope and P48 light curves are present ed in the 
standard lLandoltl (|1992D and ISmith et al.l (|2002f) systems 
in Fig. [21 Note that the Swope light curves are without 
host galaxy template subtractions, and no host galaxy or 
Milky Way dust extinction corrections have been applied. 
From Fig. [TJ the site of iPTF13ebh does not appear to 
have substantial host galaxy contamination. Hence, the un¬ 
subtracted photometry should be adequate for the deter¬ 
mination of basic photometric parameters. This is corrob¬ 
orated by the comparison of the Swope and P48 r-band 
light curves. On nights when both Swope and P48 r-band 
photometry is available, the Swope points are consistently 
brighter, but by only ~ 0.01 mag on average. 

The NIR Y JH light curves were obtained using Retro¬ 
Cam, which was moved from the 1-m Swope to the 2.5-m du 
Pont telescope in 2011. The single-chip Rockwell HAWAII- 
1 HgCdTe detector and 0'/201 pixel size yield a field of view 
of 3(5 on each side. Th e images were re d uced in the stan¬ 
dard manner following IContreras et al] ( 201(1) . As in the 
optical, the NIR photometry was computed differentially 
with respect to a local sequ ence of stars. T h e loc al sequence 
was then calibrated to the iPersson et al.l (11998 ) system in 
the JH bands, with the standard fields observed during 
six photometric nights. The Y band was calibrated to a 
set of yet-to-be-published magnitudes of the Persson stan¬ 
dards. The NIR light curves are also without host galaxy 
template subtractions. The RetroCam photometry is tabu¬ 
lated in Table [3] 

We also present UV uvw 1, uvm2, and uvw2 light curves 
obtained with the U ltraviolet Optical Teles cope (UVOT; 
Rorning et al.l 12005 ) on the Swift Mission (lOehrels et al.l 
2QQ4|). The photometry follows the method outlined by 
Brown et al.l (|2012i . I2014I) . which incorporates the updated 
zero points and time -dependent sensitivity corrections of 
iBreeveld et al.1 (1201 ll) . Note that host galaxy templates are 
not yet available to remove the non-SN counts. The UVOT 
photometry is tabulated in Tabled 
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Table 2. Journal of Swope e2v photometric observations 


MJD 

B 

V 

U 


g 


r 


i 

56610.13 

18.466 

(0.024) 

18.101 

(0.021) 

20.082 (0.172) 

18.230 

(0.017) 

18.027 

(0.020) 

18.309 

(0.046) 

56611.12 

17.803 

(0.020) 

17.501 

(0.015) 



17.595 

(0.013) 

17.435 

(0.018) 

17.632 

(0.029) 

56612.14 

17.332 

(0.014) 

17.134 

(0.014) 

18.265 (0.063) 

17.193 

(0.013) 

17.056 

(0.012) 

17.166 

(0.019) 

56613.19 

16.937 

(0.019) 

16.768 

(0.013) 



16.803 

(0.011) 

16.654 

(0.011) 

16.740 

(0.013) 

56614.10 

16.604 

(0.009) 

16.482 

(0.008) 

17.247 (0.023) 

16.497 

(0.008) 

16.362 

(0.007) 

16.447 

(0.011) 

56615.13 

16.327 

(0.009) 

16.184 

(0.007) 



16.192 

(0.006) 

16.057 

(0.006) 

16.147 

(0.008) 

56617.12 

15.841 

(0.006) 

15.720 

(0.007) 



15.711 

(0.006) 

15.609 

(0.008) 

15.722 

(0.016) 

56618.12 

15.656 

(0.005) 

15.560 

(0.006) 

16.256 

(0.009) 

15.539 

(0.005) 

15.441 

(0.006) 

15.603 

(0.007) 

56619.10 

15.541 

(0.006) 

15.411 

(0.009) 

16.145 

(0.009) 

15.395 

(0.008) 

15.288 

(0.009) 

15.474 

(0.013) 

56620.11 

15.442 

(0.007) 

15.352 

(0.008) 

16.112 

(0.012) 

15.304 

(0.007) 

15.210 

(0.012) 

15.486 

(0.012) 

56621.11 

15.368 

(0.007) 

15.278 

(0.006) 

16.111 

(0.029) 

15.296 

(0.006) 

15.168 

(0.006) 

15.505 

(0.018) 

56622.10 

15.343 

(0.005) 

15.216 

(0.005) 

16.142 

(0.009) 

15.202 

(0.005) 

15.109 

(0.006) 

15.423 

(0.007) 

56623.09 

15.343 

(0.006) 

15.181 

(0.008) 

16.191 

(0.011) 

15.189 

(0.005) 

15.099 

(0.009) 

15.469 

(0.011) 

56624.07 

15.361 

(0.007) 

15.167 

(0.006) 

16.270 

(0.013) 

15.185 

(0.005) 

15.068 

(0.006) 

15.480 

(0.010) 

56625.09 

15.401 

(0.006) 

15.176 

(0.006) 

16.397 

(0.014) 

15.211 

(0.005) 

15.059 

(0.007) 

15.507 

(0.009) 

56630.08 

15.956 

(0.011) 

15.388 

(0.012) 

17.190 

(0.026) 

15.622 

(0.010) 

15.322 

(0.012) 

15.841 

(0.016) 

56631.09 

16.124 

(0.007) 

15.430 

(0.006) 

17.421 

(0.021) 

15.703 

(0.006) 

15.360 

(0.007) 

15.846 

(0.013) 

56632.09 

16.280 

(0.008) 

15.509 

(0.007) 

17.492 

(0.015) 

15.818 

(0.005) 

15.430 

(0.008) 

15.883 

(0.011) 

56633.08 

16.442 

(0.011) 

15.607 

(0.008) 

17.711 

(0.018) 

15.953 

(0.009) 

15.462 

(0.010) 

15.895 

(0.013) 

56635.08 

16.781 

(0.014) 

15.794 

(0.007) 

18.116 

(0.030) 

16.252 

(0.006) 

15.558 

(0.006) 

15.895 

(0.009) 

56636.08 

16.925 

(0.014) 

15.860 

(0.009) 

18.266 

(0.036) 

16.392 

(0.008) 

15.612 

(0.009) 

15.867 

(0.013) 

56637.07 

17.107 

(0.021) 

15.958 

(0.008) 

18.307 

(0.047) 

16.514 

(0.008) 

15.653 

(0.006) 

15.859 

(0.008) 

56638.10 

17.224 

(0.017) 

16.050 

(0.008) 

18.535 

(0.034) 

16.678 

(0.010) 

15.693 

(0.008) 

15.862 

(0.011) 

56640.06 

17.447 

(0.023) 

16.206 

(0.014) 

18.668 

(0.047) 

16.880 

(0.012) 

15.804 

(0.008) 

15.886 

(0.009) 

56641.09 

17.553 

(0.031) 

16.295 

(0.011) 

18.941 

(0.084) 

17.012 

(0.015) 

15.890 

(0.007) 

15.928 

(0.009) 

56644.05 

17.815 

(0.026) 

16.557 

(0.011) 

19.145 

(0.092) 

17.279 

(0.019) 

16.131 

(0.008) 

16.050 

(0.010) 

56645.05 

17.872 

(0.019) 

16.650 

(0.010) 

18.997 

(0.050) 

17.354 

(0.014) 

16.232 

(0.008) 

16.123 

(0.013) 

56646.05 

17.939 

(0.013) 

16.694 

(0.010) 

19.157 

(0.030) 

17.392 

(0.013) 

16.317 

(0.011) 

16.220 

(0.013) 

56647.05 

17.971 

(0.014) 

16.764 

(0.012) 

19.224 

(0.029) 

17.440 

(0.012) 

16.399 

(0.011) 

16.318 

(0.013) 

56648.05 

18.079 

(0.015) 

16.830 

(0.010) 

19.307 

(0.035) 

17.507 

(0.014) 

16.473 

(0.011) 

16.397 

(0.014) 

56649.05 

18.089 

(0.015) 

16.876 

(0.011) 

19.376 

(0.039) 

17.529 

(0.012) 

16.539 

(0.012) 

16.506 

(0.014) 

56651.07 

18.179 

(0.018) 

17.017 

(0.013) 

19.405 

(0.038) 

17.634 

(0.014) 

16.694 

(0.013) 

16.634 

(0.018) 

56652.09 

18.222 

(0.023) 

17.061 

(0.013) 

19.407 

(0.068) 

17.700 

(0.015) 

16.760 

(0.017) 

16.824 

(0.026) 

56653.04 

18.251 

(0.017) 

17.121 

(0.013) 

19.436 

(0.028) 

17.765 

(0.017) 

16.818 

(0.017) 

16.767 

(0.023) 

56654.04 

18.274 

(0.020) 

17.154 

(0.012) 

19.442 

(0.031) 

17.761 

(0.013) 

16.824 

(0.012) 

16.808 

(0.018) 

56656.04 

18.408 

(0.027) 

17.273 

(0.020) 

19.509 

(0.053) 

17.842 

(0.016) 

16.913 

(0.014) 

16.911 

(0.019) 

56657.07 

18.360 

(0.016) 

17.276 

(0.015) 

19.634 

(0.041) 

17.840 

(0.015) 

16.971 

(0.015) 

16.976 

(0.016) 

56659.06 

18.330 

(0.018) 












Notes. The MJD column lists the Modified Julian Date of each observation. The Swope photometry is without host subtraction. 


2.2. Spectroscopic observations 

Through the iPTF collaboration, low-resolution optical 
spectra were obtained from a variety of instruments. 
These include the Intermediate dispersion Spectrograph 
and Imaging System (ISIS) on the William Herschel Tele¬ 
scope (WHT), the Himalaya Faint Object Spectrograph 
Camera (HFOSC) on the Himalayan Chandra Telescope 
(HCT), the Dual Imaging Spectrograph (DIS) on the ARC 
3.5-m telescope, the cross-dispersed spectrograph FLOYDS 
on the robotic Faulkes Telescope North (FTN), the An- 
dalucia Faint Object Spectrograph and Camera (ALFOSC) 
on the Nordic Optical Telescope (NO T), the Double Beam 
Spectrograph 1DBSP: fOke fc Cunnlll982ll on the Palomar 
200-inch telescope (P200), the Deep Imaging Multi-Object 
Spectrograph (DEIMOS; iFaber et all 1200311 on the Keck 
II tele scope, the Low Resolution Imaging Spectrometer 
('LRIS: IOke et al.ll995h on the Keck I telescope, and the In- 
amori Magella n Areal Camera and Spectrograph (IMACS; 
iDressler et al.l 12011 1 on the Magellan telescope. A journal 
of observations is presented in Table [5j and the spectra are 
plotted in Fig [3] The spectra are made pu blic via the WIS- 
eREP database rtYaron fc Gal-Yaml 12012 0. 


The NIR spectrosc opic time series of iPTF13ebh was 
obtained using FIRE (Shncoe^^Mjj^OlS) on the Magel¬ 
lan telescope, GNIRS (lElias et al.1 Il998h on the Gemini 
North Telescope, and SpeX ( Ravner et al.l l2003l l on the 
NASA Infrared Telescope Facility (IRTF) as part of a joint 
CSP-CfA Supernova Group program to obtain time-series 
NIR spectroscopy of supernovae. The combination of classi¬ 
cally scheduled time on Magellan and target-of-opportunity 
queue observing on Gemini North allowed spectra to be ob¬ 
tained in a regular 2-4 day cadence when the supernova was 
young. The observing log is presented in Table [6] The NIR 
spectroscopic time series is presented in Fig. 01 

The FIRE spectra were obtained in the high-throughput 
prism mode with a 0"6 slit. This configuration yields a con¬ 
tinuous wavelength coverage from 0.8 to 2.5 pm with a res¬ 
olution of R ~ 500 in the J band. When acquiring the 
supernova, the slit was oriented along the paral lactic a ngle 
to minimize the effect of differential refraction dFilippenkol 
Il982f ). At each epoch, several frames were obtained using 
the conventional ABBA “nod-along-the-slit” technique and 
the “sampling-up-the-ramp” readout mode. The per-frame 
exposure time was between 95.1 and 158.5 seconds depend¬ 
ing on the brightness of the supernova. These exposure 
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Table 4. Journal of Swift UVOT photometric observations 


MJD 

uvw2 

uvm2 

uvwl 

56610.9 



19.805 (0.187) 

56614.0 

19.647 (0.256) 


18.515 (0.133) 

56614.9 

19.162 (0.136) 

20.773 (0.326) 

18.017 (0.101) 

56615.7 

19.277 (0.147) 

20.339 (0.250) 

17.667 (0.092) 

56617.8 

18.761 (0.116) 

19.870 (0.181) 

17.117 (0.077) 

56619.3 

18.492 (0.110) 

19.522 (0.159) 

17.016 (0.077) 

56621.6 

18.551 (0.108) 

19.465 (0.146) 

17.228 (0.079) 

56624.6 

18.708 (0.212) 


17.411 (0.083) 

56627.5 

18.835 (0.169) 

19.520 (0.225) 

17.731 (0.126) 

56629.8 

19.181 (0.185) 

19.987 (0.290) 

17.927 (0.130) 

56631.3 

19.291 (0.219) 

20.027 (0.321) 

18.138 (0.158) 

56633.3 

19.872 (0.194) 

20.214 (0.226) 

18.215 (0.110) 

56636.4 

20.009 (0.228) 

20.327 (0.249) 

18.785 (0.149) 

56642.4 

20.706 (0.341) 


19.094 (0.172) 

56648.3 

20.445 (0.277) 

20.398 (0.309) 

19.475 (0.215) 

56651.3 


20.800 (0.334) 

19.726 (0.256) 


Notes. The MJD column lists the Modified Julian Date of each 
observation. The UVOT photometry is without host subtrac¬ 
tion. 


Table 5. Journal of optical spectroscopic observations 


UT Date 

MJD 

Instrument 

^max(-S) 

^exp 

2013-11-15 

56611.00 

WHT + ISIS 

-11.9 

3.2 

2013-11-16 

56611.79 

HCT + HFOSC 

-11.1 

3.9 

2013-11-20 

56616.25 

ARC + DIS 

-6.7 

8.4 

2013-11-22 

56618.30 

FTN + FLOYDS 

-4.6 

10.4 

2013-11-23 

56619.12 

NOT + ALFOSC 

-3.8 

11.3 

2013-11-24 

56620.30 

FTN + FLOYDS 

-2.6 

12.5 

2013-11-26 

56622.41 

P200 + DBSP 

-0.5 

14.6 

2013-11-29 

56625.48 

Keck + DEIMOS 

+2.6 

17.6 

2013-12-02 

56628.32 

Keck + LRIS 

+5.4 

20.5 

2013-12-20 

56646.88 

NOT + ALFOSC 

+24.0 

39.0 

2013-12-31 

56657.00 

Magellan + IMACS 

+34.1 

49.2 


Notes. The MJD column lists the Modified Julian Date of each 
observation. We adopt the time of B maximum of JD 2,456,623.4 
and the time of explosion of JD 2,456,608.4 for the phase relative 
to B-band maximum, f m ax(B), and phase relative to explosion, 
texp, respectively. 


Table 6. Journal of NIR spectroscopic observations 


UT Date 

MJD 

Instrument 


texp 

Tint 

2013-11-14 

56610.14 

FIRE 

-12.8 d 

2.3 d 

52.8 

2013-11-16 

56612.17 

FIRE 

-10.7 d 

4.3 d 

21.1 

2013-11-20 

56616.10 

FIRE 

-6.8 d 

8.2 d 

25.4 

2013-11-23 

56619.23 

GNIRS 

-3.7 d 

11.4 d 

40.0 

2013-11-26 

56621.96 

GNIRS 

-1.0 d 

14.1 d 

32.0 

2013-11-30 

56626.08 

FIRE 

+3.2 d 

18.2 d 

19.0 

2013-12-04 

56630.43 

GNIRS 

+7.4 d 

22.7 d 

30.0 

2013-12-09 

56635.06 

FIRE 

+12.1 d 

27.2 d 

25.4 

2013-12-14 

56640.05 

FIRE 

+ 17.1 d 

32.2 d 

16.9 

2013-12-20 

56646.05 

FIRE 

+23.1 d 

38.2 d 

16.9 

2013-12-24 

56650.32 

SpeX 

+27.4 d 

42.5 d 

30.0 

2013-12-27 

56653.03 

FIRE 

+30.1 d 

45.2 d 

16.9 

2014-01-01 

56658.03 

FIRE 

+35.1 d 

50.2 d 

16.9 


Notes. The MJD column lists the Modified Julian Date of each 
observation. We adopt the time of B maximum of JD 2,456,623.4 
and the time of explosion of JD 2,456,608.4 for the phase relative 
to B-band maximum, f ma x(B), and phase relative to explosion, 
texp, respectively. T n t represents the total on-target integration 
time in minutes. 


UT date 

11 Nov 21 Nov 01 Dec 11 Dec 21 Dec 31 Dec 



rest-frame phase relative to B-band maximum (days) 

Fig. 2. UV, optical and NIR light curves of iPTF13ebh. The 
uvw2, uvm2 and uvw 1 light curves are obtained with UVOT on 
Swift; the uBVgri light curves are obtained with the e2v imager 
on Swope; and the Y JH light curves are obtained with Retro- 
Cam on du Pont. All of these are without host subtractions. For 
the r-band light curve, open red square symbols represent host- 
subtracted r-band data from P48. A downward arrow marks the 
date of a non-detection imag e from P48 down to a magnitude 
limit of 21.7 mag. SNooPy fits (Burns et al.ll20lH 'l are plotted for 
the more sparsely-sampled NIR light curves. Blue and red verti¬ 
cal dash lines mark the dates when the optical and NIR spectra 
are taken, respectively. Two black vertical dotted lines are drawn 
at the inferred time of explosion and time of B-band maximum. 
The absolute magnitude is computed using a dist ance modulus 
of 33.63, derived from the host recession velocity of Woods et al . 

: (l 2006 i) 


times were chosen such that an adequate signal was ob¬ 
tained in each frame without saturating the bright sky lines 
in the K band. At each epoch, an AOV star was observed 
close to the science observations in time, angular distance 
and air m ass for tellur i c corr ection, as per the method de¬ 
scribed in lVacca et al.1 (120031) . 

The data were reduced using the IDL pipeline 
firehose, specifically designed for the reduction of FIRE 
data. The pipeline performed steps of flat fielding, wave¬ 
length calibration, sky subtraction, spectral tracing and ex¬ 
traction. Th e sky flu x was modeled using off-source pixels as 
described bv lKelsoii ( 2003 ) and subtracted from each frame. 
The spectral ext raction was then performed using the op¬ 
timal technique ( Hornefll986 :). a weighting scheme that de¬ 
livers the maximum signal-to-noise ratio while preserving 
spectrophotometric accuracy. Individual spectra were then 
combined with sigma clipping to reject spurious pixels. Cor- 
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rest wavelength (pm) 

0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 



observed wavelength (pm) 


Fig. 3. Optical spectra of iPTF13ebh. The UT date of obser¬ 
vation, phase relative to explosion and phase relative to B-band 
maximum are labeled for each spectrum. The gray vertical bands 
mark the regions of the strongest telluric absorptions. 

rections for telluric absorption w ere performed using the 
IDL tool xtellcor developed by IVacca et al.1 (12003 ). To 
construct a telluric correction spectrum free of stellar ab¬ 
sorption features, a model spectrum of Vega was used to 
match and remove the hydrogen lines of the Paschen and 
Brackett series from the A0V telluric standard. The result¬ 
ing telluric correction spectrum was also used for flux cali¬ 
bration. 

The GNIRS spectra were observed in the cross-dispersed 
mode, in combination with the short-wavelength camera, a 
32 lines per mm grating, and 0V675 slit. This configuration 
allows for a wide continuous wavelength coverage from 0.8 
to 2.5 pm, divided over six orders and yields a resolution 
of R ~ 1000. The observing setup was similar to that de¬ 
scribed for FIRE observations. Because of the higher resolu¬ 
tion for GNIRS, higher per-frame exposure times, between 
240 and 300 seconds, were chosen. The slit was positioned 
at the parallactic angle at the beginning of each observa¬ 
tion. An A0V star was also observed after each set of science 
observations for telluric and flux calibration. 

The GNIRS data were calibrated and reduced using the 
XDGNIRS pipeline, specifically developed for the reduction of 
GNIRS cross-dispersed data. The pipeline is partially based 
on the REDCAN pipeline for reduction of mid-IR imaging and 
spectrosco py from CANARICA M on t he Gran Telescopio 
Canarias (jGonzalez-Martm et al.l 12013 ). The steps began 
with pattern noise cleaning, non-linearity correction, locat¬ 
ing the spectral orders and flat-fielding. Sky subtractions 
were performed for each AB pair closest in time, then the 


rest wavelength (pm) 

0.8 0.9 1.0 1.1 1.2 1.4 1.6 1.8 2.0 2.2 2.5 



observed wavelength (pm) 


Fig. 4. NIR spectra of iPTF13ebh. The UT date of observation, 
phase relative to explosion and phase relative to B-band max¬ 
imum are labeled for each spectrum. The gray vertical bands 
mark the regions of the strongest telluric absorptions. 


2D spectra were stacked. Spatial distortion correction and 
wavelength calibrations were applied before the ID spec¬ 
trum was extracted. To perform telluric correction, the stel¬ 
lar hydrogen lines were first removed from the telluric star 
spectrum. The IRAF0task telluric then interactively ad¬ 
justed the relative wavelength shift and flux scale to di¬ 
vide out telluric features present in the science spectrum. 
A blackbody spectrum was then assumed for the telluric 
star for the flux calibration. Finally, the six orders were 
joined to form a single continuous spectrum. 

The SpeX spectrum was obtained in the single-prism 
low-resolution (LRS) mode. In this configuration, SpeX 
yields an uninterrupted wavelength coverage from ~ 0.7 
to 2.5 pm and a resolution of R ~ 200. The same ABBA 
“nod-along-the-slit” technique was used to obtain 12 frames 
of 150-second exposures, totaling 30 minutes of on-target 
integration time. The details in the observation and re duc¬ 
tion procedures can be found in lGushing et al.l (|2004f) and 
iMarion et al.l (120091) . 


2 The Image REduction and Analysis Facility (IRAF) is dis¬ 
tributed by the National Optical Astronomy Observatories, 
which is operated by the Association of Universities of Research 
in Astronomy, Inc., under cooperative agreement with the Na¬ 
tional Science Foundation. 
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3. Photometric properties 

The photometric properties of a SN la serve as useful in¬ 
dicators of possible peculiarities. We use the preliminary 
Swope light curves (without host galaxy subtractions) and 
pre- and post-discovery images from P48 to establish some 
basic light curve parameters for iPTF13ebh. 

The earliest data points in the r-band light curve from 
P48 and Swope allow for an accurate explosion time esti¬ 
mation. We use a Markov Chain Monte Carlo code to deter¬ 
mine the time of explosion and the associated uncertainty. 
The exponent is treated as a free parame ter to al low for de¬ 
viations from the t 2 fireball regime fe.g.. IZhener et al1l20l4 
iGoobar et al.l 120151) . The time of explosion for iPTF13ebh 
is determined to be between the earliest detection and the 
non-detection image the night before, on 2013 November 
11.85 UT or JD eX pi os i on = 2,456,608.4 ± 0.2. That is 0.59 
day after the non-detection, 0.36 day before the first detec¬ 
tion, 1.30 days before the discovery image, and 2.29 days 
before the first NIR spectrum. The result is robust against 
varying the number of light curve points included in the cal¬ 
culations. The resulting exponent is 1.5 ±0.1, substantially 
different from the fireball regime and the average exponents 
determined using large samples of SN l a light curves (e. g., 
IConlev et al.ll2006l : lHavden et al.l[201~0l : iFirth et al.ll2014 ). 

In order to test the effect of host subtraction on the 
time of explosion, we subtract the 180-degree rotated im¬ 
age of the host from each Swope r-band image. The host, 
NGC 89 0 , is morp hologically classified as a SAB galaxy 
(Ide Vaucouleurs et al.1 fl99ll) . and its projected shape in 
the plane of the sky is assumed to have axial symmetry 
(Fig. [[]). The resulting subtracted r-band P48 and Swope 
light curve is again analyzed with the Markov Chain Monte 
Carlo code. The time of explosion, 2,456,608.2 ± 0.2, is 
fully consistent with what we found using the subtracted 
P48/un-subtracted Swope light curve. This points to the 
importance of the earliest two P48 points for the determi¬ 
nation of time of explosion. 

As a comparison to measuring the explosion time with 
the r-band light curve, we also constrained the explosion 
time with a fit of the v oc f~ 0 22 power law of lPiro fe Nakaij 
(l2013h to the Sill A0.6355 /zm velocity time evolution 
(Fig. [5]). If this velocity-inferred explosion time is earlier 
than the one measured from the light curve, it may indi¬ 
cate that the supernova had a “dark phase” where it was 
initially too dim to observe because of the delay between 
the time of explosion and the time w hen the h eating from 
56 Ni first reached the outer ejecta dPiro fc Nakarll2013l) . The 
best fit explosion time with the —0.22 exponent is 19.2±0.2 
days before R maximum. As shown in Fig. [5l the —0.22 
exponent yields an excellent fit to the data. Treating the 
exponent as a free parameter yields a consistent explosion 
date and the same exponent of —0.22. This inferred ex¬ 
plosion time implies a dark phase of ~ 4 days and that 
the 56 Ni is concentrated to the center of the ejecta. Note 
that the uncertainty in the explosion date only reflects the 
goodness of the power-law fit to the data and does not take 
into account any ap p roxi mations used to obtain Equation 
13 in lPiro fe Nakaij (j2013l) . It is however interesting to note 
that this dark phase is considerably longer than those of 
SNe 2009ig, 2011fe, and 2012cg considered in lPiro fe Nakaij 
< 111 . 

Several light-curve parameters are determined here us¬ 
ing the Swope optical and RetroCam NIR light curves. The 


7 8 9 10 15 20 25 30 



time relative to explosion (days) 

Fig. 5. The Sin A0.6355 ii m velocity time evolu tion and the 
v oc t~ 0 ' 22 power law fit of iPiro fe Nakaij (120131 1. The —0.22 
exponent yields an excellent fit to the data. The inferred explo¬ 
sion time is 19.2 ± 0.2 days before B maximum. 


R-band light curve decline rate, Amis(R) <Phillipdl 11993 1 
is a luminosity indicator, but is also found to be an indica¬ 
tor o f the strength of so me spectral fea tures in the optica l 
(e.g.- lNugent et alii 199511 and NIR (e.g.. iHsiao et al.l[2013l) . 
The color-stretch parameter sbv, also a luminosity indica¬ 
tor, is found to better d iscriminate among the fast-declining 
events than Amis (R) <Burns et al.ll2014D . Since the Swope 
light curves are so densely sampled and have good data 
quality, we elec t to i nterpolate the lig ht curves using Gaus¬ 
sian processes (IRasmussen fc Williams! l2006 i. and measure 
Ami 5 (R), SBV) an d the time of B maxim um di rectly . They 
are listed in Table [7] SNooPy fits (|Burns et al.l 1201 ill yield 
consistent values for these parameters. The time of R-band 
maximum determined her e is consistent with that found by 
IWvrzvkowski et ah (201 4). but with much higher precision. 
We also measured the host color excess and the host/Milky 
Way extinction-corrected peak magnitudes using SNooPy 
colo r-stretch fits a nd the intrinsic color loci presented in 
iBurns et al.l <2014h . These values are summarized in Ta¬ 
ble [7] 

The Atois(R) value of 1.79 places iPTF13ebh in the 
category of a fast-declining SN la, but its photometric and 
spectroscopic properties are not quite as extreme as sublu- 
minous or 91bg-like SNe la. Unlike 91bg-like objects, the 
primary maxima of the NIR iYJH light curves peak be¬ 
fore the R maximum. The secondary maxima are present 
in all iYJH bands. And its R and V absolute magni- 
tu des fall withi n the w idth-luminosity relations presented 
in iPhillips et al.1 f 1999 t). It is therefore categorized as a 
“transitional” event. In this category, iPTF13ebh is one of 
the best observed and one of the most extreme members, 
pushing toward the boundary between “transitional” and 
subluminous/91bg-like objects. We explore the subject of 
transitional events in more details and propose concrete 
defini tions for this group in Section [6] 

In iMilne et al.l ( 20131) , an apparent bimodality was ob¬ 
served based on UV-to-optical colors and UV spectroscopic 
properties. Two main groups were identified to be “NUV- 
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Table 7. Summary of basic properties of iPTF13ebh 


q (J2000) 

02 h 21 m 59?98 

S (J2000) 

+33°16'13'. , 7 

TD ° 

J -L - 7 explosion, r 

2456608.4 ± 0.2 

JDexplosion,Si n v 

2456604.2 ± 0.2 

JD max (B) 

2456623.4 ±0.1 

Amis(B) 

1.79 ±0.01 

SBV 

0.63 ±0.02 

host 

NGC 890 

heliocentric redshift 

0.0133 

distance modulus' 1 

33.63 ±0.18 

E(B - V) host 

0.05 ±0.02 

^u.max 

15.29 ±0.05 

d 

TYlB, max 

14.68 ±0.05 

WV,max 

14.62 ±0.04 

TTlg, max 

14.60 ±0.04 

TTlr, max 

14.64 ±0.04 

TTli, max 

15.11 ±0.03 

^Y,max 

14.90 ±0.02 

TTlJ, max 

14.88 ±0.02 

m H, max 

15.04 ±0.01 

M-u, max 

-18.34 ±0.19 

Mb. max. 

-18.95 ±0.19 

My,max 

-19.01 ±0.18 

•^g,max 

-19.03 ±0.18 

M-r, max 

-18.99 ±0.18 

M-i, max 

-18.52 ±0.18 

My, max 

-18.73 ±0.18 

Mj ,max 

-18.75 ±0.18 

Mh. max 

-18.59 ±0.18 


Notes. V) Derived from r-band light curve. ^ Deriv ed from the 
fit of the v oc t~ 0 ' 22 power law of lPiro fe Nakail (120131 ) to the Si II 
AO.6355 (jm velocity time evolution. V) The distance modulus 
is derived from the host recession velocity which is corrected 
for the influence of the Virgo cluster, the G reat Attractor, and 
the Shapley supercluster dMould et al.ll2000j ). The error includes 
uncertainties from peculiar velocity. The peak magnitudes 
includes reddening correction for the foreground Milky Way and 
the host extinction. 


blue” and “NUV-red”. It was also noted that the vast ma¬ 
jority of the SNe la with optical carbon detections (classi¬ 
fied based on optical CII lines) belong in the “NUV-blue” 
group. iPTF13ebh would have been classified as a “narrow- 
peaked” SN la and thus excluded from the analysis of nor¬ 
mal SNe la bv lMilne et all (|2013| ). Nonetheless, we examine 
the UV color as it appears to be a strong indicator of the 
appearance of carbon. Due to the lack of UV spectra, K- 
corrections are not appl ied to our UV ligh t curves, follow¬ 
ing the same analysis as lMilne et al.1 ( 2013 ). The uvw 1 — V 
color curve of iPTF13ebh is quite red, staying well above 
1.0 mag before B maximum (Fig. (Gj). This certainly puts 
iPTF13ebh in the “NUV-red” territory. However, note that 
the uvwl — V color curve of iPTF13ebh evolves faster than 
that of a typical “NUV-red” SN la (Fig. [B]) • This is the sub¬ 
ject of ongoing work by Milne et al. (in preparation). 

iPTF13ebh may provide a rare example of carbon detec¬ 
tion in a “NUV-red” SN la. The optical CII of iPTF13ebh is 
weak and disappears quickly (Section 14. 31) . while the NIR 
C I is strong early (Section 14.211 . It weakens quickly, but 
persists until maximum light. In a typical SN la spectro¬ 
scopic follow up in the optical, starting more than a week 
past the explosion date, the unburned carbon would have 



rest-frame phase relative to B-band maximum (days) 

Fig. 6. The umul — V color curves of iPTF13ebh, SNe 2005cf 
and 201 lby. SNe 2005cf and 201 lby are examples of the “NUV- 
red” and “NUV-blue” groups, respectively. The uvwl — V color 
of iPTF13ebh is similar to that of the “NUV-red” group, but 
evolves faster. 

been missed. Carbon studies based on optical C II features 
would then introduce a bias. Whether this bias results in 
the observed association of unburned carbon with “NUV- 
blue” SNe la requires a larger sample of NIR spectra to 
reach a firm conclusion. 

4. NIR Carbon detections and evolution 

The pre-maximum spectra of iPTF13ebh were analyzed 
using the automate d spectrum synthesis c ode SYNAPPS 
(iThomas et al.1 feoill) . derived from SYN0W dBranch et all 
120051 ). SYNAPPS uses a highly parameterized, and therefore 
fast spectrum synthesis technique, useful for identifying the 
ions that form the observed features. It has been employed 
successfully to i dentify the Cll A0.6580 p m lines in early 
optic al spectra ([Thomas et al.1120071 120111 : iParrent et all 
I2012D and the Cl A 1 .0693 /an li nes in NIR spectra taken 
near maximum light dHsiao et al.ll201.it iMarion et al.ll2015l) . 

4.1. Carbon detections in the NIR 

In the two earliest NIR spectra of iPTF13ebh (—12.8 and 
— 10.7 days with respect to B maximum), the Cl A1.0693 
/jm feature is prominent. SYNAPPS was employed to con¬ 
firm the identification and to analyze each NIR Cl line 
in detail. Several NIR Cl lines are identified at the veloc¬ 
ity of —13,000 kms -1 , representing the base of the line¬ 
forming region and are presented in Fig. [7] Cl A0.9093 pm 
is listed as the s econd strongest NIR Cl line in Table 5 
of IMarion et al.l (|2009l) and it shows up as weak notches 
in both spectra. The Cl AA0.9406, 1.1754 pm lines appear 
to be present in the spectra, but some blending makes the 
identification less secure. Cl A1.4543 pm appears isolated 
from other ions; however, the proximity to the strong tel¬ 
luric absorptions between the J and H band makes this fea¬ 
ture less than ideal for the search of unburned carbon. CI 
A1.0693 pin. on the other hand, forms one of the strongest 
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velocity (10 3 km/s) 

- 20-10 0 - 20-10 0 - 20-10 0 - 20-10 0 - 20-10 0 



rest wavelength (gm) 


Fig. 7. NIR Cl lines of iPTF13ebh in the two earliest spec¬ 
tra. The phase relative to the time of B maximum is labeled 
for each spectrum. The best-fit SYNAPPS models for both epochs 
yielded a velocity of 13,000 kms" 1 for C I, representing the base 
of the line-forming region. The best-fit as a whole and the iso¬ 
lated contributions from the main ions responsible for the spec¬ 
tral features in each wavelength region are plotted. A vertical 
band marks that velocity for each Cl line. Cl A1.4543 pm for 
iPTF13ebh is located near the region of strong telluric absorp¬ 
tion which is identified with an Earth symbol. 

features in these early NIR spectra, even comparable to the 
strength of the Mgn A1.0927 pm feature, which is usually 
the most prominent feature in this wavelength region for 

SNe la. _ _ 

Echoing past works, such as lMarion et al. (2006)), we can 
confirm here that CI A1.0693 pm is the best NIR feature for 
carbon detection. All the other CI lines shown here overlap 
with other lines or telluric features. For a second line to 
confirm the presence of C I, we recommend CI A1.1754 pm. 
This feature is detected in both iPTF13ebh and SN 1999by 
I Hoflich et al.ll2002f) . and tentatively identified in SN 20life 
( Hsiao et al ■ 20Kit). It is the second strongest CI line present 
in the NIR spectra of SN 1999by. The line blending effect 
with the On A1.1667 pm feature is expected to be small, 
as neutral oxygen is the dominant ionization state and the 
Oil A1.1667 pm line strength is intrinsically weak. 

4.2. Comparison to other SNe la 

So far, NIR CI has been detected in only four SNe la; these 
are plotted in Fig. [8] for comparis on. The first s uch d etec- 
tion is in the 91bg-like SN 1999bv. iHoflich et al.1 (l2002h pre¬ 
sented NIR spectra of SN 1999by with strong Cl A1.0693 
pm that persisted through maximum light. The morphol¬ 
ogy of the feature is similar to that of iPTF13ebh, but 
the Cl A1.0693 pm feature of iPTF13ebh weakens rapidly 
and all but disappears by maxi mum lig ht. T he sub lumi- 
nous delayed detonation model of lHofiich et al.l (l2002i! also 
predicts a very strong Cl A0.9406 pm line. For the NIR 
spectra of SN 1999by, the CI A0.9406 pm line falls outside 


velocity (10 3 km/s) 

- 20-10 0 - 20-10 0 - 20-10 0 - 20-10 0 



rest wavelength (pm) 


Fig. 8. The NIR Cl A1.0693 pm line for four SNe la with NIR 
CI detection. Each panel presents the time evolution for each 
supernova with the phase relative to B maximum labeled for 
each spectrum. The solid green curves are the best-fit SYNAPPS 
models. Dashed red and dotted blue curves show the isolated 
contributions from Cl A1.0693 pm and Mgn A1.0927 pm, re¬ 
spectively. The velocity axes on the top are plotted with respect 
to the Cl A1.0693 pm line. 


the observed wavelength region so that the prediction was 
not tested. The NIR spectra of iPTF13ebh include the Cl 
A0.9406 pm feature, but as seen in Fig. [71 the feature is not 
strong. 

SN 2011fe, discovered within one day of the explo- 
sion, has been shown to be a proto- t ypical SN la (e.g. , 
iNueent et al.ll201lllParrent et al.ll2012f) . lHsiao et al.l (|2013f) 
interpreted the flattened emission wing of Mg n A1.0092 pm 
as due to the presence of CI A1.0693 pm. It was shown to in¬ 
crease in strength toward maximum light (Fig. [5]). We note 
that the profiles and the behaviors of the CI features of 
SN 2011fe and iPTF13ebh are drastically different. While 
the CI line strength increases toward maximum light for 
SN 2011fe, the Cl line in iPTF13ebh is strong, but de¬ 
creases in strength rapidly soon after explosion. iPTF13ebh 
is the only one of the four SNe la with detected Cl that 
shows this behavior. We note that while the first two spec¬ 
tra of iPTF13ebh exhibit Cl profiles similar to those of 
SN 1999by, subsequent spectra show only traces of Cl, 
manifested in the flattened wing of Mg n A1.0092 pm, much 
like the Cl line profile of SN 2011fe. 

SN 2014J developed a very similar Cl li ne pro fi le to 
that o f SN 2011fe. Indeed, as mentioned by iHsiao et al.l 
( 2013 1). the flattened wing of Mg II A1. 0092 p m i s a common 
characteristic of normal-bright SNe Ia. lMarion et al.l (|2015l) 
showed a possible detection of CI in the SYNAPPS fit of 
the maximum-light NIR spectrum. Here we take advantage 
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of the densely-cadenced observations of SN 2014J to show 
the evolution of the C I feature. The SYNAPPS fits in Fig. [5] 
suggest a moderate increas e in strength. Unfortunately, the 
late d iscovery ofSN 2014J dZheng et al.ll20 1 4 iGoohar et al.l 
l2014f) did not permit early spectroscopic follow up that is 
required to show definitiv ely the increase i n CI strength, as 
was shown in SN 2011fe (iHsiao et al.ll2013 ). 

SN 2014J has been shown to be a normal SN la, with 
signif icant reddening associated with a large amount of 

dust (|Amanullah et al1l20l4lFolev et aO20l4lBrown et al.l 

I2014T). and re lativel y high velocity for Sin AO.6355 pm 
dMarion et alJ 1201511 that places the supernova near the 
boundaries b etween the normal and high-velocity sub¬ 
classes (e.g., iBenetti et alJ 120051: IW ang et ail [20091. In 
large samples of optical spectr a, iFolatelli et al.l (1201211 and 
ISilverman fe Filippenkol (l2012t) both observed that SNe la 
with carbon tend to have lower Sill velocities, while the 
objects without carbon span the entire range of Si II veloc¬ 
ities. This phenomenon may be a result of the difficulty in 
detecting Cn A0.6580 pm at high velocities where it over¬ 
laps with the Sin A0.6355 pm line. The detection of NIR 
Cl A1.0693 pm in the moderately-high velocity SN 2014J 
shows that searching for signatures of unburned carbon in 
the NIR may remedy this observational bias. We explore 
the differences in NIR Cl and optical Cn lines further in 
the following subsection. 

In Fig. [9j we compare the CI velocity evolution for the 
four SNe la with CI detections. As the location of the CI 
line minimum in most cases cannot be directly measured 
(Fig. [5]), we plot the Cl velocities, representing the base of 
the line-forming region from the best-fit SYNAPPS models. 
They therefore have no associated uncertainties. For both 
iPTF13ebh and SN 1999by, the CI velocity evolution is flat 
(although, SN 1999by does not have early-phase data to 
rule out an early decline). The line-forming region probes 
the inner carbon-rich layers. For SNe 2011fe and 2014J, the 
C I velociti es cont inu e to decline until only a few days before 
maximum. iParrent et al.l (120111) presented the C il-to-Si II 
velocity ratio measured from a sample of optical spectra 
and showed flat time evolution and a ratio of ~ 1.1 for the 
majority of the SNe la studied. iPTF13ebh certainly falls 
outside of this norm with varying C 1-to-Si n velocity ratio 
increasing from ~ 1 to ~ 1.3. 

4.3. Comparison of NIR Cl and optical Cn 

The search for unburned material in SNe la has largely 
been focused at optical wavelengths and principally on the 
CII A0.6580 pm feature, which disappears quickly a few 
days past explosion. The recent detections of the NIR CI 
A1.0693 pm line, still present near maximum light, suggest 
that going to the NIR is a less biased way of investigating 
unburned material. In this section, we compare the NIR CI 
and optical CII features for the four SNe la with reported 
NIR CI detection. 

In Fig. [lOl we compare the NIR C I and optical C II 
features of iPTF13ebh as inferred by the SYNAPPS mod¬ 
els. As mentioned in the previous subsection, the NIR C I 
of iPTF13ebh is the strongest ever observed, but weakens 
rapidly in dramatic contrast to the behavior of the NIR C I 
lines of SNe 2011fe and 2014J. The first two NIR spectra, 
taken at —12.8 and —10.7 days relative to B maximum, 
both show a strong Cl A1.0693 /nn line, and there is evi¬ 
dence for a weak C n A0.6580 pm line in the first two optical 
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Fig. 9. The time evolution of the Cl A1.0693 pm velocities for 
the four supernovae with CI detections. Various symbols repre¬ 
sent different objects as noted. Blue, black, and red symbols rep¬ 
resent SNe la with optical light-curve decline rates in the range 
Ami5(B)<1.0, 1.0<Ami5(B)<1.6, Amis(.B)>1.6, respectively. 
The velocities are taken from the Cl velocity representing the 
base of the line-forming region in the best-fit SYNAPPS models. 
The Sill A0.6355 pm velocity of iPTF13ebh is also plotted for 
comparison. 
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spectra taken at —11.9 and —11.1 d ays relative to B rnax i- 
mum. In the classification scheme of lFolatelli et alJ (12012D . 
the Cn A0.6580 pm feature would be classified as “F” for 
flat Sill AO.6355 pm emission, and not as “A” for a defini¬ 
tive carbon absorption. The optical CII line also weakens 
rapidly. Both the weak optical C n and the strong NIR CI 
features appear at approximately the same epoch. It is pos¬ 
sible that the optical C n feature was stronger at an even 
earlier epoch, then later yielded to the growing Cl due to 
recombination. 

We now consider the optical CII lines for other SNe la 
with NIR CI detections. The optic al spectra of SN 1 999by 
from iGarnavich et al.l (120041) and iMatheson et al.l (|2008l) 
showed no strong CII features from —5 da ys throu gh to 
maximum light, while the NIR spectra of iHoflich et al.l 
(2002) showed strong detections of several CI features dur¬ 
ing the same phases. The opti cal CII A0.6580 pm fea¬ 
ture was detected in SN 201 lfe (iParrent et al.l 12012 4. and 
i s weaker th an th e NIR Cl, especially at maximum light 
dHsiao et al.ll20T! . The NIR Cl is weak early on but grad¬ 
ually increases in strength up until maximum light. While 
Cl was de tected in several p re-maximum NIR spectra of 
SN 2014J dMarion et alJl2015h . CII was not detected in the 
optic al spectra taken during the same phases (iGoobar et al.l 
I2014D . These optical CII and NIR CI comparisons are sum¬ 
marized in Fig. m along with the SYNAPPS fits and the 
illustration of the CI and CII contributions. 

While the line profile, strength and time evolution of 
the NIR Cl lines are drastically different in iPTF13ebh, 
SN 1999by, and normal SNe la like SNe 2011fe and 2014J, 
the common theme is: except at the very early epochs, 
the strength of the NIR CI lines is always observed to be 
much stronger than that of their optical CII counterparts 
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Fig. 10. Comparison between the optical CII A0.6580 /im and the NIR CI A1.0693 /im lines of the four SNe la with CI detections. 
The optical/NIR pair is selected such that the spectra are close in phase for each SN la. The phases relative to B maximum are 
labeled. The velocity axes are plotted with respect to each carbon line. Since the two earliest optical spectra of iPTF13ebh are 
taken less than one day apart and have identical Sill/Cn line profiles, they are plotted together. The SYNAPPS fit was done on the 
day —11.1 spectrum. Except at the very early epoch of SN 2011fe, the NIR C I line is always stronger than the optical C II line. 


(Fig. [TOD. The laboratory line strength of C II A0.6580 pm is 
a few orders of magnitude stronger than that of CI A1.0693 
at a reasonable range of temperature. The observations 
of stronger NIR CI lines indicate that SN la, regardless of 
their luminosity, produce conditions in the outer layers that 
are favorable to a much higher abundance of CI than CII. 
This appears t o be in contradiction with the prediction by 
iTanaka et al.l (l2008f) . that the dominant carbon ionization 
state in the outer parts of the ejecta is Cll. 

Four Cl detections is obviously too small a sample 
to reach a firm conclusion. However, we note here that 
the small number of CI detections is not due to the lack 
of or weak Cl in SNe la, but the the small sample size 
of pre-maximum NIR spectra. The most recent analyses 
of SN la NIR spec troscopy have all r evealed Cl A 1.0693 
p m to be presen t (iHsiao et al.l 1201 3: iMarion et all 12015 : 
IStritzinger et al.ll2015D 

5. Other Spectroscopic Properties 

Besides probing unburned material in the ejecta, early 
time-series NIR spectra provide mea sures of several phys ¬ 
ical para met ers i n SN la explosions (I Wheeler et al .111 99811 . 
IHsiao et, al.l (1 201 .‘ill introduced quantitative measurements 
for the IT-band break and Mg II velocity. We explore these 
quantities and other spectroscopic properties of iPTF13ebh 
in this section. 

5.1. Optical Spectroscopic Properties 

Pseudo-equivalent widths (pEW) have been widely adopted 
as a tool for quantifying supernova spectral features. We 
measure the pEW of several optical spectral features of 
iPTF13eb h nea r ma ximum light, in the same fashion as 
iFolatelli et alJ (120131) . In particular, the pEW measure¬ 
ments of the Sin AA0.5972, 0.6355 pm lines are 48.9 ± 0.6 
A and 125.2 ± 0.5 A, respectively. First introduced by 


iNugent et alJ ( 1995f) . the I?(Siii) ratio is defined as the 
ratio of the depth of the Sill AA0.5972, 0.6355 pm lines, 
and was shown to correlate with the peak luminosities of 
SNe la. A similar correlation had been shown between the 
pEW ratio of the same features and light- curve decline rate 
(|Blondin et alJl2012t IFolatelli et afll2013l) . The Sin pEW 
ratio of iPTF13ebh fits well with other fast-declining SNe la 
in this correlation. 

Using the definition given by IFolatelli et alJ (120131) . 
iPTF13eb h is plac ed firmly in the “cool” category of 
iBranch et all (2006;), for pEW{Si II 0.5972) > 30 A. The 
maximum light absorption complex near 0.43 pm is at¬ 
tributed to Mg II and Fe III in normal SNe la. This feature 
is dominated by strong Ti I I ab sorpti ons in the spectra of 
91bg-like events. IFolatelli et alJ (120131) named the pEW of 
this feature pEW2> and defined “extreme cool” (eCL) events 
to have pEW3 > 220 A. This group of SNe la is largely 
comprised of 91bg-like events, where heavy Tin absorp¬ 
tions increase pEW3 drastically. For iPTF13ebh, pEW 3 
is measured to be 104 ± 1 A, confirming the weakness of 
Ti II. SYNAPPS fits are also performed on the pre-maximum 
optical spectra of iPTF13ebh. The complex near 0.43 pm 
described above is attributed to Mg II, Fe III, and to a lesser 
extent, Sim (Fig, fill) . No Till contribution is required to 
account for the observed profiles. 

High-velocity Ca II triplet absorption features have been 
a topic of interest. Interpreted as a density enhancement in 
the outer region, the presence of high-velocity Can could 
be a signature of circumstellar material bound in the pro¬ 
genitor system which was swept up during the explosion 
(iGerardv et al.ll2004l) . iMazzali et al. ( 200 5) observed that 
almost all SNe la with early spectra show high-velocity 
Call. The question is: for the few SNe la that do not show 
high-velocity Ca II, would they show high-velocity Ca II if an 
early spectrum were available? iPTF13ebh is discovered ex¬ 
ceptionally young with the first optical spectrum obtained 
3.2 days past explosion. (The Call triplet is present in the 
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Fig. 11. SYNAPPS fit of the early optical spectrum of iPTF13ebh, 
taken at 11.1 days before B maximum. Isolated contributions 
from ions included in the fit are also shown. 


NIR spectra, but the detached high-velocity component, 
if present, would be outside of the wavelength coverage.) 
SYNAPPS fits for all pre-maximum optical spectra show that 
only a photospheric component is required to yield excel¬ 
lent fits to the observed Call triplet near 0.8 pm. (See 
Fig. HU for example.) For iPTF13ebh, there is no evidence 
of a detached high-velocity Call feature as early as 3.2 days 
past explosion. This is in agreement with previous results 
that low-luminosity SNe la are less likely to show strong 
detached high-velocity Ca l l features at early ti mes (|Hsiaol 
120091 : IChildress et al.lf20l4 iMaeuire et al.ll2014l) . 


5.2. Delayed Detonation Model 


IHoflich e t alJ (12 0021) constructed a set of delayed detonation 
fe.g.. lKhokhlovll99lD models that span the range of normal- 
bright to subluminous SNe la. For iPTF13ebh, we select the 
model with a deflagration-to-detonation transition density 
of p tr = 16 x 10 6 g cm -3 , which synthesized 0.268 M 0 of 
56 Ni. The choice is based on the peak absolute brightness 
and decline rate of iPTF13ebh, with no further fine tuning 
to match the model to the observed spectra. The model 
yields a peak absolute magnitude in B of —18.22 mag, a 
decline rate in B of Ami 5 (i?)= 1.82, and the rise time 
between the time of explosion and B maximum of 15.4 days, 
comparable to the observed values for iPTF13ebh (Table[7]). 

For the calculations of the explosions, light curves, and 
spectra, we use the radiation transport code HYDRA, 
which includes hydrodynamic solvers, Eddington tensor 
and Monte Carlo methods for low-energy photon, gamma- 
ray and positron trans port, and tim e -dependent, nuclear 


and atomic networks 


s port, an 
(iHoflichl 


1990 ; Hoflich et alJ 119951 : 


lHoflichir2002l 120091 1 Penney fe Hoflich l2014f) . The synthetic 


light curves and spectra were re com puted using updated 
atomic data from De_et_aL d2010|) and Friesen et al. (I2014T) 
as des cribed in iDiamond et alJ (|2014 ) and iTelesco et al l 
(l2015h . 

In Fig. H21 the iPTF13ebh and the pt r = 16 x 10 6 g 
cm -3 model spectra are compared between 4 and 38 days 


past ex plosion. The lin e ident ificat ions are simil a r to t hose 
given in I Wheeler et al.l (119981) and IHoflich et alJ (120021) . At 
early epochs, the spectra are dominated by blends of Cl, 
O I, Mg II, Si II, Ca II, and singly-ionized elements in the iron 
group. The strong feature at ~ 0.9 pm is identified as Mgil 
AA0.9218, 0.9244 pm, Oi A0.9266 pm, and Cl AA0.9406, 
0.9658 pm. And the feature at ~ 1.05 pm is identified as CI 
A1.0693 pm and Mgil AA1.0914, 1.0951 pm (named Mgil 
A1.0927 pm in other sections). The feature at ~ 1.35 pm is 
identified as Sill A1.3650 pm. Between 8 and 11 days after 
the explosion (top panels of Fig. fT21) . the feature at 1.05 pm 
becomes too broad with the emission wing too strong in the 
model compared to the observations. This feature is pro¬ 
duced by Mgil, with contributions from CI on the blue side 
and the suppression of the emission component due to iron- 
group elements on the red side. Increasing the 56 Ni abun¬ 
dance in the outer ejecta would have two desired effects: 
higher ionization which would lead to decreased contribu¬ 
tion from C I, and a stronger suppression on the emission 
wing from the iron-group elements. These may be achieved 
through a brighter model, stronger outward mixing of 56 Ni, 
or higher progenitor metallicity. Note that the ionization of 
carbon is sensitive to changes in transition density only for 
models in this range of brightness. 

In general, the influence of iron-peak elements increases 
with time, as the spectrum-formation region probes deeper 
ejecta layers. However, the 1.05 pm region remains rela¬ 
tively clear of iron-peak lines before maximum, demon¬ 
strating again that the C I and Mg II lines here are ideal 
for studying unburned material and the boundary between 
carbon and oxygen burning. At later time, around and after 
B maximum (bottom panels of Fig. fT?!) . the model features 
have much lower velocities compared to iPTF13ebh, indi¬ 
cating that the model has a line-forming region that recedes 
too fast. This suggests a need for an increase in the model 
56 Ni production. Alternatively, the outward mixing of 56 Ni 
may increase the opacity of layers of partial burning. Note 
that using the the explosion date inferred from a power-law 
fit to the Sin A0.6355 pm velocity measurements improves 
the situation, but the evolution of the NIR spectral features 
is still much too fast in the model spectra. 

At early time, less than a week past explosion, the model 
produces CI lines that are slightly weaker than ones ob¬ 
served in iPTF13ebh. The ionization balance of C I and C II 
depends sensitively on the recombination rate and heating. 
CI could become stronger by enhanced density as a result 
of interaction wit h circumstellar material in the progen¬ 
itor system (ICerardv et al.l l2004h which increases the re- 
combinat i on rate, o r by cooling due to early CO formation 
(IHoflich et al .11 19951 ). The strength of the model C I lines can 
also be increased by lowering the 56 Ni mass. However, the 
fast-evolving model spectra mentioned above suggests that 
lowering the model 56 Ni would cause further discrepancy 
between the model and observed spectra. 

5.3. H-band Iron-peak Feature 

The 17-band iron-peak feature, a complex formed by 
Fell/Coll/Nin, is the most prominent spect ral feature for 
a norm al SN la in the NIR. First noted by iKirshner et all 
i l!)73 i. it is fortuitously located in the H band, in between 
the two strong telluric regions in the NIR, such that it is 
relatively well documented for nearby objects. As the in¬ 
terest in SN la cosmology in the NIR grows, understanding 
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Fig. 12. Comparison between the NIR. spectra of iPTF13ebh and the pt r = 16 x 10 6 g cm -3 model spectra. The phases of the 
observed spectra are noted on the plot. These are relative to the explosion date, as inferred from the early-time r-band light curve. 
The model spectra are computed at these epochs. 


this spectral feature, which shifts out of the ZZ-band for dis¬ 
tant objects, becomes especially important for if- correction 
calculations dHsiao et al.ll2013t iBoldt et al.ll2014H . The ZZ- 
band iron-peak feature also holds the promise to provide 
insight into the explosion, as the strongly-variable opacity 
provides views of ver y different depths at the same epoch 
(IWheeler et al.l [l998l) 

iHsiao et alJ (120131 ) defined a quantitative measure for 
the size of this feature. They measured the peak-to-trough 
flux ratios across the “ //-band break” nea r 1.5 /im. Here, 
we reproduce Fig. 10 of lHsiao et al.l (|2013l ) in Fig. [T3l and 
now include measurements for iPTF13ebh and of the re¬ 
cently published NIR spectra of SN 2014J (jMarion et alJ 
l2015l f . The flux ratios are measured consiste ntly across the 
entire sample. As noted in lHsiao et al.l (l2013f) . if we exclude 
the data points from the peculiar SN 1999by, the rise and 
decline of the ZZ-band break ratio is remarkably uniform. 
The high-cadence observations of SN 2014J show the de¬ 
tails of the uniform monotonic rise, although the starting 
point appears to have occurred one day earlier than for the 
other SNe la. 

The budding of the ZZ-band feature for iPTF13ebh oc¬ 
curs around 3 days past maximum light, which coincides 
with several other SNe la, including SN 201 lfe. There is 
a large range of peak ratios; however, the post-peak de¬ 
cline rate of the ZZ-band break ratio appears to be quite 
uniform. Assuming linear declines, the decline rates of 


the ZZ-band break ratio for SN 2011fe and iPTF13ebh 
are identical within the uncertainties, —0.16 ± 0.01 day -1 
and —0.16 ± 0.06 day -1 , respectively, even though the 
two supernovae have drastically different lumi nosities. The 
error- weighted mean decline rate determined bv lHsiao et al.l 
(|2013f) is —0.15 ± 0.04 day -1 . This value appears to be 
robust among SNe la with a large range of luminosi¬ 
ties. SN 201 liv, another fast-declining transitional object 
(Ami 5 (B)= 1.76 ± 0.02), shows some peculiarity in the 
time evolution of the flux ratio. The flux ratio appears to 
peak earlier than the normal 12 days past maximum in¬ 
ferred from available data, however the cadence is not dense 
enough to definitively locate the peak of the ZZ-band break 
ratio. 

IHsiao et al.l : 2() 1 ■'> ) noted the strong correlation between 
the peak ZZ-band brea k ratio and Am^AB). In Fig. [141 
we reproduce Fig. 11 of lHsiao et al.l (120131) using the same 
technique, with the addition of iPTF13ebh and a few other 
SNe la. We also plot the ZZ -band rati o a gainst the color- 
stretch parameter sbv (see iBurns et alJ (I2014H and Sec¬ 
tion [HI) . _From Fig. m and following the same procedure 
as IHsiao et~ahl ifooiijh we assume that the time evolution 
of the ZZ-band ratio reaches its peak at 12 days past B max¬ 
imum and declines linearly thereafter. For objects with two 
or more post-peak measurements, the directly-measured de¬ 
cline rate is used to extrapolate to the peak. For SNe with 
one post-peak measurement, we assume the error-weighted 
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rest-frame phase relative to B-band maximum (days) 

Fig. 13. Time evolution of the FF-band break ratio. Differ¬ 
ent symbols represent measurements of different supernovae, as 
noted in the plot. Blue, black, and red symbols represent SNe la 
with optical light-curve decline rates in the range Ami5(B)<1.0, 
1.0<Ami5(B)<1.6, Ami5(B)>1.6, respectively. The inset illus¬ 
trates the definition of the FF-band break ratio with the spectrum 
of SN 2011fe at 12 days past maximum. The linear fits to the 
post-peak decline are plotted in solid lines when more than two 
data points are available and in dashed lines when only a single 
data point is available and the error-weighted mean decline rate 
is assumed. There is a large range of peak FF-band break ratios, 
but the post-peak declines are remarkably uniform. 


mean decline rate of —0.15 ± 0.04 day -1 . The correlation 
in Fig. [TJ] remains quite strong for both light-curve param¬ 
eters, Ami5(i3) and sbv 5 with the addition of a few more 
objects. 

The correlation between the FF-band feature and light- 
curve decline rate (and luminosity) was first pointed out by 
iHsiaol (j2009t ). Using principal component analysis (PCA) 
on eight published NIR spectra, the PCA model showed 
that the SN luminosity is not only correlated with the 
emission strength of the FF-band feature, but also with 
its width. Low-luminosity and fast-declining SNe la have 
weaker FF-band features and have the iron-peak material 
responsible for the feature confined at lower velocities. 
Withi n the Chandrasekhar mass delayed d etonation sce¬ 
nario (|Khokhlovlll99ll : lYamaoka et al.lll992f) . these correla¬ 
tions are expected, as the strength and width of the FF-band 
featu r es are indicators of t he amount of 56 Ni (|Wheeler et al.l 
119981 : IHoflich et al.l [20021) . On the other hand, for explo¬ 
sion scenarios which produ ce a spread of explosion masses 
such as dynamical mergers (llben fc Tutukovll984tiWebbinkl 
ll 98411 . or in the case of strong mixing, the correlations are 
expected to be weak. 

The FF-band break ratio measurements greatly simplify 
the analysis of the complex profile shape and time evolu¬ 
tion of the FF-band feature. In Fig. [I5j we plot the time 
series of the FF-band feature of four SNe la for compar¬ 
ison. SN 2011fe shows the FF-band break forming at ap¬ 
proximately 3 days past maximum and rapidly growing in 
strength. The characteristic double-peaked profile, formed 
by Fe Il/Co Il/Ni II, is typical of a normal SN la. iPTF13ebh 



Am 15 (B) s BV 

Fig. 14. The peak of FF-band break ratio versus the optical 
light-curve decline rate Amis(B) and color stretch sbv- SNe la 
with two or more observations in the post-peak decline are plot¬ 
ted with filled circles and solid error bars, while the rest are 
plotted with open circles and dotted error bars. The correlations 
for both parameters are strong. 


and SN 2011iv are both low-luminosity, fast-declining ob¬ 
jects in the “transitional” category, yet the profile evolu¬ 
tion of their FF-band features is subtly different. At 1 — 2 
weeks past maximum, the local minimum near 1.65 pm is 
more pronounced in SN 2011iv than in iPTF13ebh. The 
lack of a local minimum near 1.65 pm is a characteris¬ 
tic of the subluminous model presented by IHoflich et al.l 
([20021 ). By ~18 days past maximum, the FF-band breaks 
for both iPTF13ebh and SN 2011iv have shifted redward 
in wavelength, due to the decreasing size of the emission 
region at late times. The latest spectra shown in Fig. [15] 
for iPTF13ebh, SNe 1999by and 2011iv share very similar 
profile shapes, with the appearance of a strong Fell/Coll 
feature near 1.5 pm, which is usually weak in normal-bright 
SNe la. SN 1999by evolved to this profile shape approxi¬ 
mately two weeks earlier than iPTF13ebh and SN 2011iv. 

Increasing central density causes a larger fraction of ma¬ 
terial to undergo electron capture during deflagration and 
produces increasing amounts stable 58 Ni (lHoflichlf2006D . Us¬ 
ing the time evolution of the detailed profile shape of the 
FF-band complex, it is possible to isolate features produced 
by stable and unstable isotopes of iron-group elements and 
extract their mass ratio. This will be a subject of future 
studies with a larger data set. 


5.4. Magnesium velocity 


Magnesium is a product of explosive carbon burning, and 
not oxygen burning. This makes magnesium lines very sen¬ 
sitive to the location of the inner edge of carbon burn¬ 
ing in velocity space (IWheeler et al. 19981 [Marion et aD 


1200 It IHoflich et alJ 120021 : Marion et al. 20091) . The strong 
and isolated NIR Mgn A1.0927 pm line is ideal for such a 
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Fig. 15. Comparison of the 5-band iron-peak feature between 
four SNe la. The phase with respect to 5-band maximum is 
labeled for each spectrum. 
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Fig. 16. The time evolution of the Mgn A1.0972 /jm veloc¬ 
ity. Various symbols represent different objects as noted on 
the plot. Blue, black, and red symbols represent SNe la with 
optical light-curve decline rates in the range Amis (5) <1.0, 
1.0<Ami5(5)<1.6, Ami5(B)>1.6, respectively. The Mgn 
A1.0972 /im velocity of iPTF13ebh decreases monotonically, in¬ 
dicating a deeper-reaching carbon burning layer than normal 
SNe la. 
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layer than normal SNe la. In this respect, iPTF13ebh is 
similar to SN 1999by. 


measurement. Wit h early NIR spectroscopy of SN 2011fe, 
iHsiao et al.1 (' 2013 : showed that the Mg II velocity declines 
rapidly and flattens between —10 to 10 days past maxi¬ 
mum. The flattened Mg II velocity appears to be ubiquitous 
among normal SNe la. This means that it is possible to use 
a single “snapshot” NIR spectrum between —10 to 10 days 
past maximum to derive a characteristic minimum velocity 
for the given object. 

We use t he same P CA model of the Mgn profile pre¬ 
sented in IHsiao et al.l (|2013l ) and the same fitting proce- 
dure using t he no nlinear least-squares fitting program of 
iMarkwardtl (2009) to determine the Mg II velocities and 
the uncertainties for iPTF13ebh and SN 2014J. These 
data points are added to Fig. 8 of IHsiao et al.l (|2013l) 
and presented in Fig. [TG] The high-cadence observations 
of SN 2014J yielded remarkably flat Mg II velocity, simi¬ 
lar to that of SN 2011fe. The Mg II velocity of SN 2014J 
is moderately high compared t o other SNe la. Not e that 
the Gaussian fitting method of dFolatelli et al.l [2013f) gives 
systematic differences in velocities of ~ 300 kms -1 . This 
illustrates the importance of using a consistent technique 
within a comparison sample. 

At early phase, the Mgn velocity of iPTF13ebh shows 
a rapid decline similar to SN 20life. In fact, because of the 
short rise time of iPTF13ebh, the rapidly declining portions 
of iPTF13ebh and SN 2011fe would line up if Fig. 17771 was 
plotted in terms of explosion date. However, while the veloc¬ 
ity evolution of SN 201 lfe flattens to reveal the base of the 
carbon burning layer, the Mg II velocity of iPTF13ebh con¬ 
tinues to decline. The Mg II line never reaches a minimum 
velocity while visible, indicating a deeper carbon burning 


6. Transitional supernovae 

The term “91bg-like,” in the strictest sense, refers to 
spectroscopically peculiar SNe la that exhibit strong 
Tin features in their maximum-light op tical spectra, 


tin teatures m tneir maximum-li gnt_oj 

like SNe 199 1bg (iFilippenko et al.l 1199^ 


_ _ _ _^ _ _, _ (Bran ch et al. 

1993 : lLeibundeut et al.lll993l) and 1999bv ( Garnavich et al 


200 4). These objects are found to be subluminous, have 


exceptionally fast-declining light curves, and are t hought 
to o riginate from old population progenitors (|Howelll 
12001 ). Subsequent studies revealed an apparent bimodal¬ 
ity among fast-declining events wh en their NIR photomet¬ 
ric properties are examined (e.g. , iKrisciunas et al.l 120091 : 
iFolatelli et alJ 120101 : iPhilliusI l2012h . Events with their pri¬ 
mary NIR maxima peaking after the R-band maximum 
are subluminous in their absolute magnitudes in all bands 
compared to normal SNe la, even after taking the width- 
luminosity relation into account. These subluminous su¬ 
pernovae also tend to lack or have very weak secondary 
NIR maxima. The terms “91bg-like” and “subluminous” are 
spectroscopic and photometric definitions, respectively. So 
far, all SNe la with late-peaking NIR maxima have been 
found to have the spectroscopic peculiarity of strong Till, 
although, not all SNe la with Tin signatures have late- 
peaking NIR maxima, as we will point out next. 

When defining “transitional” events between the fastest- 
declining normal SN la population and the 91bg- 
like/subluminous events, we chose the timing of the NIR 
primary maximum rather than the presence of Till as the 
defining feature. A “transitional” event is then defined as a 
fast-declining SN la with its NIR (iY JHK) primary max- 
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Table 8. Optical light-curve parameters and the presence of 
Ti nof notable “transitional” events 


Name 

Ami 5 (B) 

sbv 

Tin 

Data 

SN 1986G 

1.62 ±0.02 

0.65 ±0.05 

91bg-like 

(a) 

SN 2003gs 

1.79 ±0.04 

0.46 ± 0.03 

91bg-like 

(b) 

SN 2003hv 

1.48 ±0.04 

0.76 ± 0.07 

No Till 

( c ) 

SN 2004eo 

1.42 ±0.03 

0.83 ± 0.03 

No Till 

CSP; (d) 

SN 2005am 

1.53 ±0.02 

0.75 ±0.02 

No Till 

CSP; (d) 

SN 2007on 

1.89 ±0.01 

0.55 ±0.02 

No Till 

CSP; (e) 

SN 2009an 

1.64 ±0.04 

0.86 ± 0.06 

No Till 

(f) 

SN 201 liv 

1.76 ±0.02 

0.63 ±0.01 

No Till 

CSP; (g) 

SN 2012ht 

1.30 ±0.04 

0.86 ± 0.03 

No Till 

CSP 

iPTF13ebh 

1.79 ±0.01 

0.63 ± 0.02 

No Till 

CSP 


Notes. The Amis(B) and sbv values are measured us¬ 
ing the data from the quoted sources. ^ iPhillios et, al.1 
l|l987fl ^ iKrisciunas et all 1120090 C) iLeloudas et, all |2009) 
G1 IContreras et all ll201fli C) IStritzinger et, all (l201lh 
03 ISahu et al.1 llifflidll (s) Gall et al. (in preparation) 
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imum peak ing before its B-ba nd maximum. The case of 
SN 1986G dPhillips et al.l[l987h is the main driver behind 
this choice of definition. For SN 1986G, Tin is present in 
its spectra, like SN 1991bg, but its NIR absolute magni¬ 
tudes put it with the normal SN la group. It is “91bg- 
like” spectroscopically, but not “subluminous” photometri¬ 
cally. Its NIR maxima also peaks before B maximum. In 
our definition, SN 1986G is categorized as “transitional” 
with the normal SNe la and not with the “subluminous” 
group. More recently, SN 2003gs is also categorized as “tran - 
sitional” and shows strong Tin ilKrisciunas et al.112009 ). 
Readers should keep in mind that “91bg-like” and “sublumi¬ 
nous” are, strictly speaking, not the same group. However, 
since our discussions have been focused on spectroscopic 
properties, and our only NIR spectroscopic reference in the 
91bg-like and subluminous groups is SN 1999by, we use the 
two terms interchangeably. 

From Fig. [3 the primary maxima in iY JH light curves 
of iPTF13ebh clearly peak before its B maximum. And, un¬ 
like SNe 1986G and 1991bg, iPTF13ebh shows no evidence 
of Till in its optical spectra (Section 13.111 . This is clearly 
a transitional event and has one of the fastest light-curve 
decline rates in this class. For comparison, the decline rates 
of some of the more notable transitional events are listed in 
Table [U 

While the width-luminosity relation iPhillipsi 1 199311 
holds well for a large range of luminosities and Am^B), 
there is a notable break at A?7ii5(B)>1.6, where SNe la 
become faint very quickly with increased Am\§(B). Possi¬ 
bly related to this is the fact that Ami^B) fails to prop¬ 
erly classify the fastest-declining objects. The reason for 
this is that the onset of the linear decline in the .B-band 
light cu rve occurs earlier than 15 days for th e fastest de- 
cliners (|Folatelli et al.1 [20101 : iBurns et al.ll201lh . This leads 
to poorly-behaved light-curve templates and a seeming bi¬ 
modality in the photometric properties of SNe la with 
Atoi 5 (B). Use of Ami 5 (B) is th erefore pro bl ematic when 
dealing with transitional objects. I Burns et all (2014!) intro¬ 
duced a new decline-rate parameter, sbVi in an attempt to 
improve the treatment of the fast-declining objects. Using 
sbV) the fastest declining events appear less as a distinct 
population and more as a continuous tail end of normal 
SNe la. 


Fig. 17. Comparison of the Swope z-band light curve of 
iPTF13ebh to the SNooPy and SiFTO templates. The best SNooPy 
fits using Amis(B) and sbv z-band templates are plotted as dot¬ 
ted and solid curves, respectively. The best SiFTO fits with “nor¬ 
mal” and “91bg-like” templates are plotted as dash and dot-dash 
curves, respectively. 


In Fig. [TT1 the Swope z-band light curve of iPTF13ebh 
is plotted in compa rison with the “normal” and “ 91bg-like” 
fits, as described bv lGonzalez-Gaitan et al.1 (l2014f). The fits 
are d one with the light-curve fitter SiFTO ( Gonlev et all 
120081) , but with two sets of light-curve templates for com¬ 
parison: one of normal SNe la and the other constructed 
from 91bg-like objects. The fitter stretches the time axis 
to fit the iPTF13ebh data from all of the optical uBVgri 
bands. The result shows poor fits in the z band for both 
“normal” and “91bg-like” templates. iPTF13ebh has an z- 
band light curve morphology that bridges the two classes, 
with transitional properties, such as the timing of the pri¬ 
mary and secondary peaks and the size of the secondary 
peak. 

Fig. [Tf] also shows comparison SNooPy fits using 
Amis(B) and sbv templates. The fits also include all op¬ 
tical uBV gri bands. Unlike SiFTO fits, SNooPy takes into 
account the changing shape of the light curves in the redder 
filter bands for supernovae with different peak luminosities. 
The Amis(B) template fit yields an z-band light curve with 
a primary maximum that peaks too early and a secondary 
maximum that peaks too late, which is characteristic of a 
brighter supernova than the z-band light curve suggests and 
is much like the “normal” SiFTO fit described above. The 
sbv template, on the other hand, provides an excellent fit 
to the observed light curve. The effect would be even more 
pronounced in the NIR Y JH bands. We chose the z-band 
light curve for illustration, as it has much better coverage 
than the NIR bands (Fig. [2]). 

The photometric properties of iPTF13ebh place it with 
the normal SNe la, but many of its NIR spectroscopic prop¬ 
erties, such as the strong C I, the profile shape of the H- 
band feature, and the Mg n velocity evolution, are similar 
to those of 91bg-like SN 1999by. SNooPy with the sbv pa¬ 
rameterization works well for iPTF13ebh and other transi¬ 
tional objects. This suggests a continuous range of proper- 
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ties across the transition. However, the sbv parameter does 
not predict the presence of Ti II (nor does Ami^{B), see Ta¬ 
ble [SJ). It is unclear whether normal and subluminous/91bg- 
like belong in distinct groups or are objects with a common 
explosion mechanism and a continuous range of observed 
properties, but iPTF13ebh provides a connection between 
these subtypes. 

7. Conclusions 

iPTF13ebh is one of the fastest-declining and best- 
observed “transitional” objects between normal-bright and 
subluminous/91bg-like events. Overall, iPTF13ebh has 
properties that are distinct from both normal and sublu- 
rninous classes. We summarize the main findings as follows: 

1. The NIR iYJH light curves of iPTF13ebh peak before 
its B maximum, and its extinction-corrected absolute 
magnitudes obey the width-luminosity relation. 

2. Several Cl lines are detected in the early NIR spectra 
of iPTF13ebh. The CI A1.0693 /iin line at 2.3 days past 
explosion (—12.8 days relative to B maximum) is the 
strongest ever observed in a SN la. 

3. The profile of the Cl A1.0693 fin i line at first resem¬ 
bles that of SN 1999by, a subluminous/91bg-like ob¬ 
ject. However, while the strong Cl line persists until 
maximum light in SN 1999by, the CI line of iPTF13ebh 
quickly weakens and takes on a profile similar to that of 
a normal SN la. 

4. The NIR Cl A1.0693 fim line of iPTF13ebh is much 
stronger than its optical Cll AO.6580 fin i line. This 
property is general among SNe la with CI detections, 
an d contradicts the prediction of stronger C n lines by 
iTanaka et al.l (120081 ). 

5. Unlike normal SNe la, iPTF13ebh shows flat CI velocity 
evolution, indicating that the base of the unburned layer 
is already reached at the earliest phases. 

6. The 17-band break ratio of iPTF13ebh fits well in the 
peak ratio versus Atoi.R-B) relation previously found by 
iHsiao et al.l (120131) . The strong correlation is consistent 
with the Chandrasekhar mass delayed detonation sce¬ 
nario, while the correlation is expected to be weak in 
dynamical mergers. 

7. By late time, the profile shape of the 17-band complex 
resembles that of SN 1999by. The composition and den¬ 
sity of the inner core of iPTF13ebh is therefore similar 
to that of 91bg-like events. 

8. Unlike in normal SNe la, the Mgn velocity of 
iPTF13ebh monotonically decreases with time indicat¬ 
ing a deeper reaching carbon burning layer than normal 
SNe la. This behavior is similar to that of SN 1999by. 

9. There is no evidence of a detached high-velocity compo¬ 
nent for the Can NIR triplet of iPTF13ebh, as early as 
3.2 days past explosion, providing a stringent constraint 
on when this high-velocity feature could have formed, if 
it ever existed. 

10. iPTF13ebh has a substantial difference between the ex¬ 
plosion times inferred from the early-time light curve 
and the velocity evolution of the Sill AO.6355 fim line, 
implying a long dark phase of ~ 4 days. 

The data set presented here is unique, as the follow 
up campaign began shortly after explosion and spanned 
a wide wavelength coverage from the UV to the NIR. 


iPTF13ebh is also the only “transitional” SN la with ex¬ 
tensive NIR spectroscopic coverage. The photometric prop¬ 
erties of iPTF13ebh suggest that it is at the fast-declining 
end of normal SNe la. Yet, the NIR spectra reveal many 
similarities to the 91bg-like SN 1999by. The rare glimpse 
of the NIR spectroscopic properties of a transitional event 
provides a connection between normal and the peculiar 
subluminous/91bg-like events, perhaps indicating a com¬ 
mon explosion mechanism. 
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